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ABSTRACT 

o , 

■ Redshift evolution of supermassive black hole mass functions (BHMFs) is in- 
vestigated up to z ^ 1. BHMFs at intermediate redshifts are calculated in two 
ways. One way is from early- type galaxy luminosity functions (LPs); we assume an 

' Mbh — -^sph correlation at a redshift by considering a passive evolution of Lgph in 

the local relationship. The resultant BHMFs (spheroid-BHMFs) from LFs of red se- 
, quence galaxies indicates a slight decrease of number density with increasing redshift at 

Q^ ■ Mbh ^ 10 ''•^"^Mq. Since a redshift evolution in slope and zeropoint of the Mbh — -^sph 

I relation is unlikely to be capable of making such an evolution in BHMF, the evolution 

i of the sphcroid-BHMFs is perhaps due mainly to the decreasing normalization in the 

' galaxy LFs. We also derive BHMFs from LFs of morphologically selected early- type 

(— I , galaxies. The resultant BHMFs are similar to those from the red sequence galaxies, 

^ i' but show a small discrepancy at z ^ 1 corresponding to an increase of SMBH num- 

ber density by ^ 0.3 dex. We also investigate how spheroid-BHMFs are affected by 
uncertainties existing in the derivation in detail. 

The other way of deriving a BHMF is based on the continuity equation for number 
^ I density of SMBHs and LFs of active galactic nucleus (AGN). The resultant BHMFs 

(AGN-BHMFs) show no clear evolution out to z = 1 at Mbh > IO^Mq, but exhibit a 
significant decrease with redshift in the lower mass range. Interestingly, these AGN- 
k> ] BHMFs are quite different in the range of Mbh ^ lO^M© from those derived by 

■ Merloni (2004), where the fundamental plane of black hole activity is exploited. 

Comparison of the spheroid-BHMFs with the AGN-BHMFs suggests that at 
Mbh ^ lO^M©, the sphcroid-BHMFs are broadly consistent with the AGN-BHMFs 
out to z ~ 1. Although the decrease of SMBH number density with redshift suggested 
by the spheroid-BHMFs is slightly faster than that suggested by the AGN-BHMFs, we 
presume this to be due at least partly to a selection effect on the LFs of red sequence 
galaxies; the colour selection could miss spheroids with blue colours. The agreement 
between the spheroid-BHMFs and the AGN-BHMFs appears to support that most of 
the SMBHs are already hosted by massive spheroids at z 1 and they evolve without 
significant mass growth since then. 

Key words: black hole physics - galaxies: elliptical and lenticular, cD — galaxies: 
evolution. 
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1 INTRODUCTION 

Recent observations provide evidence that a mass of a super- 
massive black hole (SMBH) in a galactic nucleus is tightly 
correlated with a mass or luminosity of a spheroid compo- 
nent of its host galaxy (e.g., Magorrian et al. 1998; Mar- 
coni & Hunt 2003, MH03 hereafter). The tight correla- 
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tion suggests the presence of strong evolutionary link be- 
tween SMBH and spheroid component. Using the relation, 
an SMBH mass function (BHMF) can be derived from lo- 
cal galaxy luminosity function (LF) or velocity dispersion 
function. Meanwhile, a local BHMF can also be calculated 
from cosmological evolution of LFs of active galactic nuclei 
(AGNs) by integrating the continuity equation for number 
density of SMBHs, where mass accretion onto a SMBH is as- 
sumed to power an AGN and grow the central SMBH (e.g.. 
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Cavaliere, Morrison & Wood 1971; Small & Blandford 1992; 
Marconi et al. 2004; Shankar et al. 2004). BHMFs derived 
by this method can now be more reliable than before thanks 
to updated LPs of hard X-ray selected AGNs (Ueda et al. 
2003), which is more complete to obscured AGNs. Marconi 
et al. (2004) and Shankar et al. (2004) demonstrate that the 
local densities of SMBHs and BHMFs derived with these two 
ways agree with each other, if one adopts reasonable values 
for accretion efficiency and Eddington ratio. Furthermore, 
Marconi et al. (2004) indicate that the cosmic history of 
mass accretion rate density delineates that of star-formation 
rate density, and the ratio of the latter to the former is about 
4000, which agrees with the mass ratio of a spheroid to a 
central SMBH, again strongly suggesting the co-evolution of 
SMBHs and spheroids. 

In studying the co-evolution of SMBHs and host 
spheroids in further detail, one approach is to investigate 
redshift evolution of BHMF and correlation between black 
hole mass (Mbh) and spheroid luminosity (Laph) or mass. 
While Mbh has been measured for a substantial number 
of high redshift QSOs (e.g.. Shields et al. 2003; McLure & 
Dunlop 2004), it is technically challenging to directly mea- 
sure dormant SMBHs at cosmological distances. However, 
BHMFs at high redshifts can be computed using AGN LFs 
and the continuity equation in the same way as in the lo- 
cal universe. Also, we are now in a reasonable position to be 
able to study BHMFs at intermediate redshifts by using LFs 
of early-type galaxies. Early-type galaxy LFs have recently 
been derived out to z ~ 1 with good statistics from inten- 
sive imaging surveys such as COMBO-17 (Wolf et al. 2003; 
Bell et al. 2004b). Consequently, one can compare BHMFs 
from the galaxy LFs with those from AGN LFs, which may 
provide clues to understand the co-evolution of AGNs and 
host spheroids. 

In this paper, we will derive BHMFs at intermediate 
redshifts from galaxy LFs and AGN LFs and investigate the 
redshift evolutions. When converting galaxy LFs to BHMFs, 
a correlation between Mbh and I/gph is utilized. Although 
a correlation between Mbh and bulge effective stellar ve- 
locity dispersion (ae) is claimed to be tighter than that be- 
tween Mbh and Lsph (Ferrarese & Merritt 2000; Gebhardt 
et al. 2000), measuring ae at cosmological distance is very 
hard and we thus adopt the Mbh — isph relation in this 
work. It should also be emphasized that MH03 show the 
Mbh — Laph relation is as tight as that between Mbh and 
(Te if only SMBHs whose masses are securely determined are 
used in the analysis. In addition, MH03 suggest that the in- 
trinsic scatter of the relation in B-band is as small as in NIR 
bands. One should be aware that the derivation of a BHMF 
from galaxy LFs contains several processes which may give 
significant uncertainties to resultant BHMFs. Quantifying 
uncertainties in BHMFs when calculated from galaxy LFs 
is also aimed at in this paper. 

The layout of this paper is as follows. In the next section 
(§IUi we describe the procedure to derive a BHMF from an 
early-type galaxy LF and calculate BHMFs from early-type 
galaxy LFs up to z ^ 1. We also investigate how BHMFs 
are affected by uncertainties and unconstrained parameters 
existing in this derivation. In § |21 we derive BHMFs from 
AGN LFs and examine its redshift evolution. In § 3] we 
compare the BHMFs from early-type galaxy LFs with those 
from AGN LFs and discuss the results. Throughout this pa- 



per, we adopt the cosmological model with Ho = 70 km s 
Mpc~^, Qm — 0.3 and f^A ~ 0.7 unless otherwise stated. 



2 BHMF FROM EARLY-TYPE GALAXY LF 
2.1 Derivation of BHMF 

In this paper, we mainly use the LFs obtained by Bell et al. 
(2004b) from the COMBO-17 survey. The large survey area 
with the moderate depth brings a large number of galaxies, 
and the multi-band photometry covering from 3640 A to 
9140 A with the 17 broad- and medium-band filters allows 
one to accurately determine the photometric redshifts. Both 
of these aspects are important to derive luminosity functions 
with good accuracy and the large survey area is especially 
useful to estimate its cosmic variance. They firstly investi- 
gated the rest-frame U — V vs. My colour-magnitude di- 
agram of galaxies at a certain redshift and found the red 
sequence consistent with the colour-magnitude relation well 
established for the early-type galaxy population. They select 
galaxies on the red sequence and derive their rest-frame B- 
band LFs at redshifts from 0.25 to 1.05. BeU et al. (2004b) 
also derived the _B-band LF of the local red sequence galax- 
ies using the SDSS EDR data (Stoughton et al. 2002) by 
transforming the SDSS ugr system to the standard UBV 
system. We use the Schechter functions fitted to these LFs 
in the following analyses. A part of the survey area of the 
COMBO-17 was imaged with the HST/ACS and most of the 
galaxies on the red sequence (~ 85 %) indeed show early- 
type morphology (Bell et al. 2004a). It should be noted that 
the colour selection based on the red sequence is presumed to 
exclude blue ellipticals with star formation activity and/or 
young stellar population. In addition, small bulges in late- 
type galaxies can be missed, which perhaps results in a de- 
ficiency of the light part of a BHMF. 

In order to convert these LFs to BHMFs, firstly the 
total galaxy luminosities need to be transformed to the 
spheroid luminosities using bulge-to-total luminosity ratios 
(B/Ts). B /Ts of early-type galaxies at intermediate red- 
shifts are, however, not well constrained observationally. 
Im et al. (2002) selected morphologically early-type galax- 
ies at intermediate redshifts (most of the galaxies are at 
z = 0.2-1.0) from the HST/WFPC2 data for DEEP Groth- 
Strip Survey (DGSS) based on B /Ts (> 0.4), which are 
estimated by fitting radial surface brightness profiles with 
r^''^ spheroid and exponential disk. According to their cata- 
log, the morphologically selected early-type galaxies have a 
B/T of 0.7 on average, although there is a substantial scat- 
ter. The B/Ts do not depend on redshift or luminosity in 
their sample. Therefore we adopt B/T — 0.7 in this study, 
independently of redshift and galaxy luminosity. 

Once LFs of spheroidal components are calculated from 
the early-type galaxy LFs, they can be transformed to 
BHMFs using an Mbh — J^sph relation. In this study, we 
use those derived by MH03 for their galaxies in Group 1, 
for which measurements of Mbh and Laph are considered to 
be reliable. The relation is derived for Laph in the B, J, H, 
and K bands. Although the relation in if-band is generally 
preferred because A'-band luminosity is considered to be the 
most reliable indicator of the stellar mass of a galaxy, MH03 
show the intrinsic scatter of the relation is almost indepen- 
dent of the bands for the galaxies in Group 1. Since we will 
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mainly use the B-band LPs by Bell et al. (2004b), we adopt 
the _B-band relation described as follows; 

logMsH = (1.19 ± 0.12)(log L,ph - 10.0) + (8.18 ± 0.08). (1) 

In applying the Mbh — isph relation to an LF, the intrin- 
sic scatter of the relation needs to be considered; we adopt 
A log Mbh = 0.32 according to MHOS. 

In calculating a BHMF from a spheroid LF at a high 
redshift, an evolution of the AIbh — isph relation needs to 
be considered. In this study, we basically consider only the 
effect of passive evolution in Lsph. It should be noted that 
if all spheroid components evolve only passively (with no 
growth of SMBHs), then the BHMF at the intermediate 
redshift does not change from that at z = 0. In fact, it is not 
obvious that all spheroids evolve passively at intermediate 
redshifts. Any other evolutions in galaxy LFs and the Mbh — 
Lsph relation would cause the disagreement between BHMFs 
from galaxy LF and AGN LF, which will be discussed in §|1| 
We describe the passive luminosity evolution as Mb{z) = 
Mb[z = 0) — Qz and assume Q to be 1.4. This is evaluated 
using PEGASE Ver 2.0 (Fioc & Rocca-Volmerange 1997) 
for a stellar population with the solar metallicity formed at 
2 = 4 (the age at z = is 12 Gyr) in a single starburst with 
an e-folding time of 1 Gyr. This Q value is consistent with 
the evolution of characteristic luminosity in the COMBO-17 
LFs (z ^ 0.25). 

2.2 BHMFs from early-type galaxy LFs up to 

2 ~ 1 

Figure □ shows the BHMFs calculated using the COMBO- 
17 LFs with the prescription described above. In the upper 
panel, the B-band isph of the galaxies in Group 1 by MHOS 
is plotted against Mbh. Black solid line indicates the best- 
fitting regression line to the data (i.e., the relation ai z — 0). 
Grey solid, black dashed, and black dot-dashed line is the 
relation expected at z = 0.25,0.65 and 1.05, respectively, 
when the passive luminosity evolution is considered. In the 
lower panel, BHMFs transformed from the COMBO-17 LFs 
in the rest-frame B band at these redshifts are indicated 
with shaded regions. The envelope of each BHMF shows 
the errors of Mg and <j!>* in the Schechter function fit to 
the data (Bell et al. 2004b); the upper and lower bound 
is defined by a BHMF with the largest and smallest L*g 
and (/!>* within the fitting error, respectively. The uncertainty 
in (j>* is dominated by cosmic variance (Bell et al. 2004b) 
and we adopt the larger value of the two different estimates 
by BeU et al. (2004b). Each BHMF is indicated down to 
the black hole mass corresponding to the lowest luminosity 
among the data points in the LF. The characteristic Mbh 
corresponding to Mg at each redshift is indicated by arrow 
with the same colour as the BHMF. 

The BHMFs in Figure exhibit a redshift evolution. 
Since the characteristic Mbh does not largely change with 
redshift, this evolution is perhaps due to a density evolution 
of BHMFs corresponding to the decreasing normalization of 
the red sequence galaxy LFs (Bell et al. 2004b). We will 
discuss other possibilities in S 14.21 

Currently, the COMBO-17 survey provides one of the 
largest and most useful databases to study galaxies out to 
z ~ 1 and the derived early-type galaxy LFs are there- 
fore considered to be the most reliable so far. Neverthe- 
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Figure 1. Upper panel: The correlation between Mbh ^.nd B- 
band spheroid luminosity. Black solid line indicates the relation 
at 2 = (equation fitted to the data points by MHOS for their 
galaxies in Group 1. Grey solid, black dashed, and dotted line is 
the relation expected at 2; = 0.25,0.65 and 1.05, respectively, 
when a passive evolution of spheroid luminosity is considered. 
Lower panel: BHMFs transformed from the COMBO-17 LFs at 
redshifts of 0.25, 0.65, and 1.05 are indicated with shaded regions, 
of which widths show the errors in Mg and (p* (Bell et al. 2004b). 
The lower mass cutoff of these BHMFs corresponds to the lowest 
luminosity of the data points in the original LF. The characteristic 
Mbh corresponding to Mg at each redshift is indicated by arrow 
with the same colour as of the BHMF. 



less it is worth investigating the evolution of BHMF using 
LFs determined with other data sets, in particular to see 
whether the result is sensitive to selection criterion for early- 
type galaxy. Since early-type galaxies are selected from the 
colour-magnitude diagram in COMBO-17, we investigate 
the BHMFs derived from LFs of morphologically selected 
early-type galaxies by Im et al. (2002; the HST/WFPC2 
data for DGSS are used) and by Cross et al. (2004; the 
data were taken with the HST/ACS in the guaranteed time 
observations). In these studies, early-type galaxies are se- 
lected based on the analyses of radial surface brightness 
profiles. In Im et al. (2002), the rest-frame _B-band LFs 
are derived in the two redshift bins: 0.05 < z < 0.6 and 
0.6 < z < 1.2, while in Cross et al. (2004), they are ob- 
tained at 0.5 < z < 0.75 and 0.75 < z < 1.0. 

In Figures Hand 01 the BHMFs converted from the LFs 
of morphologically selected early-type galaxies by Im et al. 
(2002) and Cross et al. (2004), respectively, are presented 
by shaded regions, showing the uncertainties due to the fit- 
ting errors of Mg and <j>* in the LFs. Again these BHMFs 
are calculated by assuming B/T — 0.7 and using the B- 
band Mbh — isph relation by MHOS corrected for passive 
luminosity evolutions. The low mass cutoff of the BHMF 
corresponds to the lowest luminosity among the data points 
in the original LF. The BHMF from the COMBO-17 LF at a 
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Figure 2. The BHMFs converted from the early-type galaxy LFs 
in the DGSS (Im et al. 2002) are indicated with shaded regions, 
whose widths represent the uncertainties of the BHMFs due to 
the fitting errors of and 0* in the LFs. The pair of solid lines 
describes the BHMF from COMBO-17 at similar redshifts and 
the separation of the two lines indicates the uncertainty of the 
BHMF. The low mass cutoff of the BHMF corresponds to the 
lowest luminosity among the data points in the LF. 



similar redshift is also indicated with solid lines, showing the 
upper and lower bounds defined by considering the errors of 
A'lg and 0* in the LF. These comparisons demonstrate that 
the BHMFs from the morphologically selected early-type 
galaxy LFs axe consistent with those from the COMBO-17 
survey. It should be mentioned that the BHMFs obtained 
from COMBO-17 LFs tend to lie below those from the mor- 
phologically selected early-type galaxy LFs at the higher 
redshifts; the discrepancy in number density of SMBHs at a 
given AIbh is estimated to be ~ 0.3 dex. This may indicate 
that the contribution of blue spheroids which are not in- 
cluded in the COMBO-17 LFs becomes larger towards z ^ 1. 

2.3 Source of uncertainty in the derivation of a 
BHMF 

In the previous subsections, we showed the procedure to 
convert galaxy LFs to BHMFs at intermediate redshifts and 
presented the resultant BHMFs up to 2 ~ 1. However, there 
are several sources of uncertainty in the derivation. Some of 
them are related to the Mbh — isph relation such as fitting 
error of the relation to the data points, uncertainty of its 
intrinsic scatter, and choice of Mbh — isph relations. Others 
are due to the fact that B /T value and passive luminosity 
evolution are not well constrained. In what follows, we will 
investigate how BHMFs are affected by these uncertainties. 
Unless otherwise noted, we will begin all the calculations to 
derive BHMFs with the LF of red sequence galaxies at 2 = 
(Bell et al. 2004b, the SDSS EDR data are used). We caU 
this input LF "iLF" hereafter. 
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Figure 3. Same as Figure|3 but for the BHMFs converted from 
the early-type galaxy LFs computed by Cross et al. (2004) using 
the HST/ACS data. 



2. 3. 1 Fitting error and uncertainty of intrinsic scatter in 
A/bh — Jjsph relation 

In Figure 0] the fitting error in the B-band AfsH — isph re- 
lation estimated by MHOS is demonstrated in the top panel. 
Three relations are shown here: One is the best-fitting re- 
gression line (solid line), and the other two are those which 
result in the most massive or least massive BHMFs within 
the ±1(7 fitting error. The BHMF using either of the three 
relations is indicated in the middle panel. In calculating the 
BHMFs, the intrinsic scatter A log Mbh = 0.32 (MH03) is 
taken into account. This shows that the part of a BHMF at 
A/bh ^ Iff'^MQ is affected by this uncertainty. 

We also calculate BHMFs for several values of intrinsic 
scatter of the relation and show the results in the bottom 
panel; (grey solid line), 0.2 (dotted line), 0.3 (black solid 
line), and 0.4 (dashed line). This affects again the massive 
end of a BHMF (Mbh > 10® '^ M©). It is suggested by MH03 
to be between 0.3 and 0.4 (see also McLure & Dunlop 2002), 
but we note that if observational errors in Mbh are under- 
estimated, the intrinsic scatter would be smaller. 

2.3.2 Mbh — ^sph relations in B hand and K band 

We will compare BHMFs derived from iLF using the A^bh — 
Lsph relations calculated by MH03 either in the B band 
or K band. In order to apply the A"-band relation to iLF 
{B band), we need to convert the isT-band relation to a B- 
band relation by correcting it for B — K colour. This can be 
performed by using the average colour of early-type galaxy 
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Marconi & Hunt (2003) 
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Figure 4. Top panel: The Mbhs and LgphS of the galaxies in 
Group 1 by MHOS are plotted and the best-fitting regression line 
to the data points is indicated by solid line. Dashed (dotted) line 
shows the relation which results in a BHMF most (least) biased 
towards the massive end within the 1 cr fitting errors, respectively, 
when applied to an LF. Middle panel: BHMFs are computed from 
iLF with the A/bh — L^ph relations demonstrated above and they 
arc indicated with the same line styles as in the top panel. The 
intrinsic scatter of the relation (AlogMen = 0.32; MHOS) is 
considered to derive the BHMFs. Bottom panel: BHMFs are cal- 
culated for several values of intrinsic scatter of the relation; 
(grey solid line), 0.2 (dotted line), O.S (black solid line), and 0.4 
(dashed line). 



population (e.g., Marconi et al. 2004; McLure & Dunlop 
2004; Shankar et al. 2004). We adopt B-K = 3.75 estimated 
by Girardi et al. (2003) for local elliptical and SO galaxies 
in field and group environments. This colour is similar to 
that used in Marconi et al. (2004), where Mz — K = 4.1 
for ellipticals and 3.95 for SOs {Mz is Zwicky magnitude) 
are used in applying the isT-band Mbh — ^sph relation to 
the LFs from the CfA survey (Marzke et al. 1994). These 
colours are based on actual measurements by Kochanek et 
al. (2001) and they are converted to B — K = 3.59 and 3.44, 
respectively, by using the equation B = Mz — 0.51 given by 
AUer & Richstone (2002). 
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Figure 5. Upper panel: B-band M-qh — L^-py^ relations. Black solid 
line is the B-band Mbh — ^sph relation by MHOS and the other 
lines indicate those obtained by correcting the i^-band relation for 
A-CMR (grey solid line), T-CMR (dotted fine), and the average 
B — K colour (dashed line) . Lower panel: BHMFs calculated from 
iLF using the above B-band relations are indicated. 



However, this correction ignores the colour-magnitude 
relation (CMR). In fact. Bell et al. (2004b) suggest that the 
red sequence galaxies in the COMBO-17 survey are on the 
CMR which is consistent with those found for E/SOs both in 
clusters (Bower, Lucey, & Ellis 1992; Terlevich, Caldwell & 
Bower 2001) and fields (Schweizer & Seitzer 1992). There- 
fore, we attempt to take this into account. Since CMRs in 
B — K have rarely been investigated, we model a CMR using 
a population synthesis code (Kodama & Arirnoto 1997) so 
as to reproduce an observed CMR in a certain set of filters 
and we then derive a CMR m. B — K and Mb using the 
model. In this modelling calculation, we follow the recipe by 
Kodama et al. (1998), where they modeled the CMR [V-K 
and Mv) measured in the Coma cluster (Bower et al. 1992) 
based on the galactic wind scenario. 

It needs to be mentioned that, while a CMR is normally 
defined using colours within a metric aperture, the one us- 
ing total colours suits better for this study. Total colours in 
more luminous (larger) galaxies are presumed to be system- 
atically bluer than those measured within an aperture due 
to more severe effects of the colour gradients (e.g., Peletier 
et al. 1990; Tamura & Ohta 2003); a given aperture can 
sample only the reddest part of a luminous (large) galaxy, 
but it can include the total light of a faint (small) galaxy. 
The slope of a CMR using total colours (T-CMR hereafter) 
is therefore expected to be flatter than that using colours 
within an aperture (A-CMR hereafter). This aperture eff'cct 
on the CMR in the Coma cluster (Bower et al. 1992) has 
been investigated by Kodama et al. (1998) and the slope of 
the T-CMR at My ^ -20 mag is estimated to be ~ 30% 
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flatter. Hence we consider a CMR with a slope flatter by 
30% than the A-CMR by Bower et al. (1992) as a T-CMR.^ 
In Figure|3 the B-band A/bh — J^sph relations are shown 
in the upper panel. Black solid line is the B-band Mbh — iaph 
relation by MHOS and the other lines indicate those obtained 
by correcting the i^-band relation for the average B — K 
colour, A-CMR, and T-CMR. In the lower panel, BHMFs 
calculated from iLF using these B-band relations are in- 
dicated with the same line styles as above. It is indicated 
that the _ft'-band relation is not fully transformed to the 
B-band relation by MHOS even if a CMR is taken into ac- 
count, and the BHMFs come towards the less massive end 
than the case where the B-band relation by MHOS is ap- 
plied to iLF. This may imply that, while the tightness of 
the B-band Mbh — isph relation is nearly the same as that 
of the AT-band relation, the relations are not equivalent to 
each other. One possible reason for this discrepancy may be 
that B band luminosity is less good indicator of stellar mass 
due to effects of dust extinction and/or young stellar pop- 
ulation for some of the galaxies used in the analysis of the 
Mbh — -t'sph relations. It should be noted that the BHMF 
comes towards the less massive end when the average B — K 
colour is used than when a CMR is considered because the 
adopted average colour is bluer than colours of luminous 
early-type galaxies. Consequently, their B-band luminosi- 
ties are overestimated and therefore smaller values of A/bh 
are assigned to spheroids with a certain B-band luminosity. 
For the same reason, the BHMF also depends on choice of 
A-CMR or T-CMR, but the difference turns out to be very 
small. 



2.3.3 Choice of Mw 
authors 
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relations from different 



Next we will compare BHMFs derived using the A^bh — isph 
relations by different authors: Ferrarese & Merritt (2000, 
FMOO), Merritt & Ferrarese (2001, MFOl), McLure & Dun- 
lop (2002, MD02)2, and MHOS. For FMOO, we adopt the re- 
lation for their Sample A. In MFOl and MD02, the relation 
is obtained in the V band and R band, respectively, and they 
are converted to B-band relations using T-CMRs modeled 
in the same way as explained earlier. Also, the zeropoints of 
the relations by MFOl and MD02 are shifted by the amounts 
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Figure 6. Same as Figure|S] but AfBH—^sph relations by different 
authors and BHMFs derived using them are compared. Upper 
panel: Black solid line, grey solid line, dotted line, and dashed 
line indicates the _B-band A^bh ~ ^sph 

relation by MHOS, FMOO, 
MFOl, and MD02, respectively. In MD02 and MFOl, the relation 
is originally derived in the V and R band, respectively, and they 
are converted to the _B-band relations using T-CMR (see text for 
details). Lower panel: BHMFs derived by applying the above B- 
band relations to iLF are indicated. The intrinsic scatter of the 
A/bh ~ ^sph relation is assumed to be 0.32 (MHOS) in all the 
calculations of the BHMFs. 



due to the differences in Bo from the value we adopt. These 
B-band relations are shown in the upper panel of Figure |S| 
with the B-band relation by MHOS. They are applied to iLF 
and the BHMFs obtained are shown in the lower panel. The 
intrinsic scatter of 0.32 around the relation is assumed in 
all the calculations. These BHMFs suggest that A/bh — L^ph 
relation depends on sample data set and consequently there 
is a substantial variation among the BHMFs. 



^ In Kodama et al. (1998), Hq = 50 km Mpc~^ is assumed 
and hence the slope of a T-CMR would be less flat in our cos- 
mology {Ho = 70 km s~^ Mpc~^). But here we aim at seeing 
how BHMFs vary by considering the CMRs and precise deter- 
mination of the slope is beyond our scope. We note that colours 
of ellipticals less luminous than My = — 20 mag are considered 
to be robust to the aperture correction and the actual T-CMR 
therefore has a break at My ~ — 20 mag (see Figure S in Kodama 
et al. 1998). We ignore this and determine the zeropoint of the 
T-CMR so as to reproduce the total colours of the ellipticals with 
My ^ —20 mag. Although this indicates that the T-CMR gives 
too red colours to less luminous ellipticals, the impact of this on a 
BHMF is very small and does not affect the following discussions. 

Strictly speaking, since about half of the sample consists of 
QSOs at 0.1 < z < 0.5 in MD02, the Mbh - Lsph relation is not 
allowed to be used here because some evolutionary effects may 
already be incorporated. In practice, however, the AfBH — i'sph 
relation derived using only the local inactive galaxies in MD02 is 
the same as that from the whole sample. 



2.3.4 Bulge-to-total luminosity ratio 

B /T at intermediate redshift is only loosely constrained 
from observations at the moment and we have to await for 
future works to see the validity of the current assumption 
(B/r = 0.7 for all the red sequence galaxies). It is therefore 
worth demonstrating BHMFs for a range of B/T to keep 
it in mind as uncertainty. In Figure Q BHMFs for B /T = 
0.5, 0.7, and 0.9 are compared. These BHMFs are obtained 
from iLF and the B-band A/bh — ibuigc relation by MHOS. 
This indicates that a BHMF is affected by choice of B/Ts 
at AfBH ^ lO*Af0. 

2. 3. 5 Model of passive luminosity evolution 

We consider a passive luminosity evolution of an old stellar 
population which formed at z = 4 with a starburst of which 
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Figure 8. Top panel: The Mbh — ^sph relations at z = 0.55 
are calculated considering passive evolution of Lgph and are com- 
pared with the relation at 2; = by MHOS (black solid line). We 
consider three models for passive luminosity evolution of old stel- 
lar population; zj- = 10 (dashed line), 4 (grey solid line) and 3 
(dotted line). Bottom panel: BHMFs at 2 = and 0.55 derived 
using the above relations are indicated with the same line styles. 
Note that all the BHMFs (including that at 2 = 0) are calculated 
from the COMBO-17 LF at 2 = 0.55 (not iLF). 

e-folding time is 1 Gyr. In this case, the luminosity evolution 
in the B band is described as Mb{z) = Mb{z = 0) — Qz 
and Q = 1.4. The Q value can be estimated using other 
models of stellar populations and its dependency on choice of 
parameters may need to be treated as uncertainty of BHMF 
at high redshift. Since formation redshift (zf) of a stellar 
population is probably the most important parameter upon 
which the Q value largely depends, we consider two other 
cases: Zf = 10 and 3. The Q value is estimated to be 1.2 and 
1.6, respectively. Using these Q values including 1.4 for Zf — 
4, we derive passively evolved B-band Mbh — isph relations 



at z = 0.55 and apply these relations to the COMBO-17 
LF at this redshift (not iLF). The results are displayed in 
Figure |H| We also calculate BHMFs at z — 1.05 using the 
LF and the passively evolved Mbh — isph relation at this 
redshift and show the results in Figure |5] 

2. 3. 6 Summary of the uncertainties 

Table Q shows a summary of uncertainties in a BHMF in- 
vestigated earlier. In each column, the possible range of log- 
arithmic SMBH number density {A\og(f>M) caused by each 
source of uncertainty is indicated at the five BH masses: 
log Mbh ~ 7.5, 8.0, 8.5, 9.0 and 9.5. The error sources in- 
vestigated are arranged in the columns (1) — (6) as follows: 
(1) the la fitting error of the B-band Mbh — isph rela- 
tion (see the middle panel of Figure 3J, (2) the uncertainty 
of intrinsic scatter (0.3 — 0.4) in the Mbh — isph relation 
(the lower panel of Figure (3) choice of filter bands {B 
or K) of Mbh — isph relation (Figure |SJ , (4) choice of the 
Mbh — Lsph relations by different authors (Figure [SJ (5) 
B /T (0.5 — 0.9; Figure 0, and (6) choice of passive evolu- 
tion models (2:/ = 3 — 10; A log (pM is calculated for BHMFs 
at z = 1.05 (Figurel^J). In calculating Alog<^M for (3), we 
compare a BHMF computed with the B-band Mbh ~ isph 
relation obtained by correcting the /('-band relation for the 
average B — K colour (i.e., CMR is not considered) with a 
BHMF derived with the B-band relation by MH03. In (4), 
the largest and smallest SMBH number densities are taken 
from the BHMFs in Figure|S| In addition to these, A log (j>M 
due to the errors of characteristic luminosity and normal- 
ization in a COMBO-17 LF is shown for comparison in the 
column (7). The errors in the LF at z = 0.65 are adopted 
to calculate A log (pM here, but the errors in a BHMF are 
similar if LFs at other redshifts are used (see the lower panel 
of Figure . 
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log A/bh 



A log 0A/ 





(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


7.5 


0.08 


0.02 


0.02 


0.11 


0.02 


0.02 


0.16 


8.0 


0.04 


0.02 


0.07 


0.26 


0.09 


0.04 


0.20 


8.5 


0.10 


0.01 


0.27 


0.37 


0.27 


0.15 


0.28 


9.0 


0.37 


0.14 


0.57 


0.50 


0.53 


0.32 


0.41 


9.5 


0.80 


0.45 


0.99 


0.96 


0.88 


0.54 


0.58 



dt 



dL 



,(3) 



Table 1. Summary of uncertainties in a BHMF at the five BH 
masses. Each number shows the possible range of logarithmic 
SMBH mass density {A\og<f>M) at a given Mbh caused by each 
source of uncertainty: (1) the Icr fitting error of the B-band 
A^BH~^sph relation (see the middle panel of Figure|l}, (2) the un- 
certainty of intrinsic scatter (0.3 — 0.4) in the A^bh — ^sph relation 
(the bottom panel of Figure|^, (3) choice of filter bands {B or K) 
of Mbh ~ ^aph relation (Figure |2J, (4) choice of the Mbh ~ ^sph 
relations by different authors (Figure [SJ, (5) B/T (0.5 - 0.9; 
Figure 171. (6) choice of passive evolution models {zj = 3 — 10; 
Alog(/)M is calculated a.t z = 1.05 ("Figure 151^. and (7) errors of 
characteristic luminosity and normalization in a COMBO-17 LF 
(the errors of the LF at 2 = 0.65 are considered here; see the lower 
panel of Figure^. In calculating Alog(f>M for (3), we compare a 
BHMF computed with the B-band Mbh ~ ^sph relation obtained 
by correcting the X-band relation for the average B — K colour 
(i.e., CMR is not considered) with a BHMF derived with the B- 
band relation by MH03. In (4), the largest and smallest SMBH 
number densities are taken at each Mbh. 



This table demonstrates that most of the uncertainties 
are negligible at Mbh ^ 10* M©, while they are significant in 
the range of Mbh 1O*M0. The uncertainty related to the 
photometric band selection and colour correction and that 
related to selection of the Mbh — isph relations by different 
studies are the most serious and they amount to an order 
of magnitude at the massive end. Those due to the fitting 
error of the S-band Mbh — isph relation and the possible 
uncertainty in B/T also appear to be substantial. 



3 BHMF FROM AGN LF 

In this section, we investigate the cosmological evolution of 
BHMFs derived ft-om AGN LPs (AGN-BHMFs hereafter). 
The results will be compared with those from COMBO-17 
LFs (spheroid-BHMFs hereafter) in § |1| By assuming that 
only mass accretion grows the central SMBH at a galactic 
centre and galaxy mergers are not important in its growth 
history, the time evolution of a BHMF 4>M{MBii,t) can be 
described by the continuity equation: 



!^M(MBH,t) 
dt 



+ 



d 



-['/'M(MBH,i){M(MBH,t)>] = 0, (2) 



where (M(Mbh, i)) represents the mean mass accretion rate 
at a given SMBH mass Mbh and at a cosmic time t. Fur- 
thermore, if we assume a constant radiative efficiency e 
(= L/M(?, where L and M is the bolometric luminosity 
and the mass accretion rate, respectively) and a constant 
Eddington ratio A (= L/L^dd, where L^dd is the Eddington 
luminosity) for all the AGNs, the second term of the above 
equation can be simply related to a (bolometric) luminosity 
function of AGNs. This finally reduces the continuity equa- 
tion to: 



L = AAfBHcVtE 



where tEdd is the Eddington time and ip{L, t) is the AGN LF 
(for details, see Marconi et al. 2004). We note that xp{L,t) 
is the number of AGNs per dlogL, while <?!)a/(Mbh, i) is 
the number of SMBHs per dMnn- Hence, once the form 
of an AGN LF is known as a function of redshift, one can 
integrate this equation to obtain BHMFs at any redshifts 
starting from the initial condition, either in time decreasing 
order (from a high redshift to z = 0) or the inverse (from 
z = to higher redshifts). 

Following the procedure adopted by Marconi et al. 
(2004), here we derive AGN-BHMFs at intermediate red- 
shifts starting from a BHMF at z = 3 as the initial condition 
(it is assumed that all the SMBHs at z = 3 were shining as 
AGNs). In the calculation, we use the hard X-ray AGN LF 
(HXLF) by Ueda et al. (2003, U03 hereafter), which is de- 
scribed by a luminosity-dependent density evolution model 
(LDDE model; see their § 5.2 for details). To take into ac- 
count the contribution of "Compton-thick" AGNs to the 
total mass accretion rate, we multiply a correction factor of 
1.6 independently of the AGN luminosity. The luminosity- 
dependent bolometric correction described in Marconi et al. 
(2004) is adopted. The Eddington ratio and the radiative ef- 
ficiency are assumed to be constant, A = 1.0 and e = 0.1, re- 
spectively, also based on the study by Marconi et al. (2004) . 
The results are shown in Figure 1101 indicating that while 
the BHMFs in the range of Mbh ^ 10** M© hardly change 
out to z ~ 1, they exhibit a clear redshift evolution at 
Mbh ^ 10** M0. That is, almost aU SMBHs with a mass 
larger than 10® M© formed at 2 > 1, while lighter SMBHs 
grow later, suggesting a downsizing of SMBH evolution. 

It needs to be pointed out that there are several uncer- 
tainties in the AGN-BHMFs thus far derived as follows: 

(i) There are ranges of values in Eddington ratio and ra- 
diative efficiency which give a reasonable fit of an AGN- 
BHMF to the spheroid-BHMF at z = 0. The distribution 
studied by Marconi et al. (2004) suggests a possible range 
of A = 0.1 - 2.0 and e = 0.04 - 0.15 within 1 a uncer- 
tainty. In Figure we exemphfy AGN-BHMFs at 2 = 
and 0.65 for five sets of A and e; either of A or e is fixed to 
the adopted value (A — 1.0 or e = 0.1) and the largest or 
smallest value within the uncertainty is chosen for the other 
parameter. This plot indicates that, as expected from equa- 
tion the normalization of an AGN-BHMF is altered and 
the AGN-BHMF is shifted along the Mbh axis by changing 
A, while only the normalization is affected by changing e. In 
the following analysis, AGN-BHMFs will be calculated for a 
number of pairs of A and e on the 1 a contour provided by 
Marconi et al. (2004; see their Figure 7) and the envelope 
of these AGN-BHMFs that gives a possible range of SMBH 
density at a given Mbh will be adopted as uncertainty of 
AGN-BHMF. 

(ii) The assumption of constant Eddington ratio and ra- 
diative efficiency for all the AGNs is perhaps too simple. Al- 
though Marconi et al. (2004) claim that the local spheroid- 
BHMF can be well reproduced by the AGN-BHMF with a 
constant A ~ 1.0 and e ~ 0.1, the solution only gives a suffi- 
cient condition to the limited constraints at z = 0. Further- 
more, Heckman et al. (2004) claim that AGNs have various 
Eddington ratios and the ratios seem to depend on A/bh 
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Figure 10. AGN-BHMF at a redshift of 0, 0.4, 0.8, 1.2 and 2.0 
is plotted with open circles, solid circles, triangles, squares, and 
crosses, respectively. These AGN-BHMFs are calculated with the 
continuity equation and the HXLFs by U03 (see text for details). 
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Figure 11. AGN-BHMFs at 2 = (upper panel) and those at 
z = 0.65 (lower panel) for five sets of X and e as shown in the 
upper panel. 



Kawaguchi ot al. (2004) also propose that supcr-Eddington 
accretion is essential for a major growth of SMBH. 

(iii) There is ambiguity in the continuity equation itself; 
if a merging process should be added as a source term in the 
equation, a resultant BHMF would change. 

(iv) Although the HXLF by U03 accounts for all the 
Compton-thin AGNs including obscured AGNs, the uncer- 
tainties in the estimate of Compton-thick AGNs directly af- 
fect the resulting BHMFs. 

Quantifying all of these uncertainties but (1) requires 
substantial theoretical works and/or new observational data 
and is beyond the scope of this paper, but one must keep 
them in mind. 




6 7 8 9 10 

Log Mg„/Mg 

Figure 12. AGN-BHMFs at a redshift of 0.1, 0.6, and 1.2 ob- 
tained in this study by integrating the continuity equation start- 
ing from 2 = 3 are plotted with solid circles in the top, middle 
and bottom panel, respectively. The error bars represent the up- 
per and lower bounds of SMBH density at a given Mbh that are 
calculated from AGN-BHMFs for the values of A and e within 
the 1 cr uncertainty. Asterisks show AGN-BHMFs obtained by 
Merloni (2004). 



It may be intriguing to compare these AGN-BHMFs 
with those recently obtained by Merloni (2004) , where a new 
method is employed to investigate the redshift evolution. 
They introduce a conditional luminosity function (CLF), 
which is the number of active black holes per unit comoving 
volume per unit logarithm of radio (Lr) and X-ray {Lx) 
luminosity and is defined so that by integrating a CLF over 
the range of Ljj or Lx, the radio LF (RLF) or X-ray LF 
(XLF) of AGNs is obtained, respectively. In order to relate 
Lr and I/x of AGN to Mbh without any assumptions on ac- 
cretion rate, they use an empirical relation among Lr, Lx, 
and Mbh (fundamental plane of black hole activity; Mer- 
loni, Heinz & Di Matteo 2003) . A CLF also needs to satisfy 
a constraint that the number density of SMBHs with a cer- 
tain mass is obtained by integrating a CLF over the ranges 
of Lr and Lx which are determined by the fundamental 
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plane. Consequently, once a BHMF is obtained at a redshift 
of z, a CLF can be computed using RLF, XLF, and BHMF 
as constraints.^ Note that RLFs and XLFs from observa- 
tions are available up to high redshifts; the RLFs obtained 
by Willott et al. (2001) and the HXLF by U03 are used in 
Merloni (2004). Given a functional form of accretion rate'*, 
a mean accretion rate can also be calculated as a function 
of Mbh- Using a BHMF and a mean accretion rate at z, a 
BHMF a,t z + dz can be derived from the continuity equa- 
tion. Likewise, BHMFs at higher redshifts can successively 
be calculated and hence RLF, XLF, and BHMF at z — 
are firstly needed. In Merloni (2004), the local BHMF is ob- 
tained using LFs of galaxies with different morphologies by 
Marzke et al. (1994) and the empirical relationships among 
spheroid luminosity, velocity dispersion, and Mbh- 

Figure ^] shows the comparison of AGN-BHMFs at 
redshifts of 0.1, 0.6, and 1.2 calculated in the two differ- 
ent methods: (I) BHMFs derived by integrating the conti- 
nuity equation from z = 3 to lower redshifts using the U03 
HXLF with A = 1.0 and e — 0.1. The error bars represent 
the upper and lower bounds of SMBH density at a given 
Mbh that is calculated from AGN-BHMFs for a number of 
pairs of A and e on the 1 a contour provided by Marconi 
et al. (2004). (II) Those obtained by Merloni (2004), which 
are integrated from 2; = to higher redshifts using the CLF 
with the fundamental- plane relation. This comparison indi- 
cates that while the agreement of the AGN-BHMFs is good 
in the massive end, the discrepancy at Mbh ^ 1O^'^M0 
is substantial at all redshifts. There could be several rea- 
sons for this disagreement. One possibility may be related 
to the choice of an initial spheroid-BHMF; we note that the 
BHMF at z = 0.1 in Merloni (2004) has a normalization ~ 
10 times smaller at Mbh ~ IO'^Mq than the local BHMF 
independently estimated by Marconi et al. (2004, see their 
Figure 2). Another reason could be the different estimate 
for the mean mass accretion rate. In fact, we find that the 
second term of the equation Q averaged between z = 0.9 
and 2; = 0.1 calculated in method I is significantly larger 
than that in method II in the range of A/bh ^ lO'^'^M©. We 
just point out the facts in this paper and leave further dis- 
cussions for future studies. In the next section, we adopt the 
AGN-BHMFs calculated by method I for comparison with 
the spheroid-BHMFs. 



4 DISCUSSIONS 

4.1 Comparison of Spheroid-BHMFs with 
AGN-BHMFs 

In Figure 1131 spheroid-BHMFs are compared with AGN- 
BHMFs up to 2: ~ 1. The spheroid-BHMFs transformed 
from the COMBO-17 LFs at redshifts of 0., 0.25, 0.45, 0.65, 
0.85, and 1.05 are indicated with shaded regions (the BHMF 
at z = was derived from the LF of the red sequence galax- 
ies in the SDSS EDR; see Appendix in BeU et al. 2004b). 



In this calculation, the fundamental plane of black holes is as- 
sumed to be independent of redshift. 

4 In Merloni (2004), -Lx/^Edd = f(M,m) where m = 
£accMc^/LE(jd is adopted (sacc is accretion efficiency). 



The widths of the shaded regions are determined by consid- 
ering not only the errors of Mg and (/>* but also the uncer- 
tainty associated with band transformation between B and 
K of the Mbh — ^sph relation, which is the most significant 
uncertainty among those investigated (see S 12. 311 . The upper 
bound of the shaded region is the BHMF derived by applying 
the B-band Mbh — isph relation by MH03 to an LF with the 
largest values of characteristic luminosity and normalization 
within the errors. To determine the lower bound, a B-band 
Mbh — isph relation calculated by correcting the A"-band 
relation for the average B — K colour of early-type galaxies 
is applied to an LF with the smallest values of characteristic 
luminosity and normalization within the errors. The trans- 
formation of the Jf-band relation to B band is performed 
only at 2: = and the B-band relations at high redshifts are 
then obtained by correcting it for passive luminosity evolu- 
tion in the B band. The spheroid-BHMF is depicted with 
shaded region down to the mass corresponding to the lowest 
luminosity among the data points in the LF. At the masses 
lower than this cutoff, the upper- and lower-bound BHMFs 
are indicated by dashed lines. The characteristic A/bh corre- 
sponding to L*g at each redshift is indicated by arrow. Two 
arrows are shown in each panel; the one at the more mas- 
sive end shows the characteristic Mbh of the upper-bound 
BHMF, and the other is of the lower-bound BHMF. The 
solid curve, which goes through the middle of the upper- 
and lower-bound BHMFs at 2 = 0.25, indicates the BHMF 
at 2 = 0.25^ calculated with a _B-band Mbh — isph rela- 
tion converted from the K-hanA relation by correcting for 
T-CMR. This BHMF is plotted in all the panels as a fidu- 
cial of comparison. The AGN-BHMFs calculated with the 
continuity equation and the HXLFs by U03 at the same 
redshifts as the spheroid-BHMFs are overplotted with open 
circles. The error bars indicate the upper and lower bounds 
of SMBH density allowing for the 1 a uncertainty of A and 
e. 

From Figure^! it is suggested that at A/bh > 10* Mq, 
the spheroid-BHMFs^ are broadly consistent with the AGN- 
BHMFs out to 2 ~ 1. This agreement between the spheroid- 
BHMFs and the AGN-BHMFs appears to support that most 
of the SMBHs are hosted by massive spheroids already at 
2^1 and they evolve without significant mass growth since 
then. The discrepancy at Mbh ^ lO^'^M© between the 
spheroid-BHMFs and the AGN-BHMFs is presumed to be 
due at least partly to the fact that small bulges in late-type 
galaxies tend to be excluded in selecting the red sequence 
galaxies and thus their contributions are not expected to be 
included in the COMBO-17 LFs or the spheroid-BHMFs. 
In fact, galaxy LFs are not well constrained down to such 
low luminosities and future observations therefore need to 



5 We adopt the BHMF not at z = but at z = 0.25 just for 
consistency; the LFs at z ^ 0.25 are calculated with the data 
from the COMBO-17 survey, while the LF at z = is from the 
SDSS data. 

^ In Figure lTi^ the spheroid-BHMF at 2: = 0.25 appears to exceed 
that at z = 0. This is because of the large increase of character- 
istic luminosity in the COMBO-17 LF from z = to 0.25. This 
increase is significantly larger than that predicted for passive evo- 
lution, while the rate of luminosity evolution at z ^ 0.25 is fully 
consistent with passive evolution out to z ~ 1. The origin of the 
large luminosity increase at the low redshift is currently unknown. 
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Figure 13. The sphcroid-BHMFs transformed from the COMBO-17 LPs at redshifts of 0, 0.25, 0.45, 0.65, 0.85, and 1.05 are indicated 
with shaded regions, whose widths are determined by considering not only the errors of and <j)* but also the uncertainties associated 
with band transformation in the derivation of a BHMF (see text for details). The spheroid-BHMF at z = was derived from the LP of 
the red sequence galaxies in the SDSS EDR (see Appendix of Bell et al. 2004b). The low mass cutoff of the shaded region corresponds 
to the lowest luminosity among the data points in the original LP. In the lower mass range than the cutoff, the upper- and lower-bound 
BHMPs are indicated by dashed lines. The characteristic Mbh corresponding to Mg at each redshift is indicated by arrow. Two arrows 
are shown in each panel; the one at the more massive end shows the characteristic Mbh of the upper-bound BHMP, and the other is 
of the lower-bound BHMP. Solid curve indicates a BHMP at z = 0.25 calculated with a B-band Mbh ~ ^sph relation converted from 
the X-band relation by correcting for T-CMR. This BHMP is displayed in all the panels as a fiducial of comparison. The AGN-BHMPs 
calculated with the continuity equation and the HXLPs by U03 at the same redshifts as those of the spheroid-BHMFs are overplotted 
with open circles. The error bars indicate the upper and lower bounds of SMBH density allowing for the 1 a uncertainty of A and e as 
in Pigure VPTl 
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be awaited for any further discussions on the discrepancies 
in the light end of BHMF. It is interesting to point out that 
while the AGN-BHMFs at Mbh ^ 10* M© do not signif- 
icantly evolve out to ^ ~ 1, the spheroid-BHMFs exhibit 
a slight redshift evolution (see also Figure and Figure 
llOjl . We note that the uncertainties of the spheroid-BHMFs 
considered in Figure 1131 are larger than those in Figure 
where the evolution in spheroid-BHMF may look clearer. 
One possible reason for the difference in evolution between 
the spheroid-BHMFs and the AGN-BHMFs is a selection 
effect on the red sequence galaxy LFs; if there are more blue 
spheroids with on-going star formation and/or young stellar 
population towards z = 1 then their contribution is more 
likely to be missed from the red sequence galaxy LF and 
the spheroid-BHMF at higher redshifts. In fact, as men- 
tioned earlier, the BHMFs obtained from morphologically 
selected early-type galaxy LFs tend to exceed those from 
the COMBO-17 LFs at 2 ~ 1 (Figures 13 and EJ . 

4.2 Does the correlation between A/bh and host 
spheroid mass evolve with redshift? 

An alternative interpretation of the possible difference in 
evolution between the spheroid-BHMFs and the AGN- 
BHMFs may be a difference at high redshift between the 
actual Mbh — ^sph relation and our assumption. In other 
words, there may be an evolution of the Mbh ~ isph rela- 
tion other than the passive luminosity evolution of spheroid. 
In order to see whether this can be the case, it is worth 
examining how a BHMF can be affected by changing the 
Mbh — isph relation. 

Here, we consider a simple case where, in the Mbh— isph 
relation (log A/bh ~ plogLsph + q), the coefficient of p or g 
varies. Figures |2l and |2| demonstrate the effect on BHMF of 
an evolution of q; a BHMF is shifted mostly in parallel to 
the Mbh axis. The effect of changing p on a BHMF is sim- 
ilar, although it can be a modification in shape of BHMF 
rather than a lateral shift. Therefore, changing neither p 
nor q moves the BHMF along the axis of the number den- 
sity of SMBH. On the other hand, the redshift evolution of 
the actual spheroid-BHMFs seems to be dominated by that 
along the number density axis, suggesting that it is diffi- 
cult to explain the slight difference in evolution between the 
spheroid-BHMF and the AGN-BHMF by changing p and q 
with redshift. 

Recently, a possible offset from the local relationship 
between Mbh and central velocity dispersion (ctq) has been 
found at z ~ 0.37 by investigating spectra of the central 
regions of galaxies hosting type 1 AGNs (Treu, Malkan & 
Blandford 2004, T04 hereafter). It is interesting to see a 
spheroid-BHMF at this redshift derived with the offset found 
by T04. In order to apply the relation at z — 0.37 to a 
COMBO-17 LF at 2: = 0.35, the Mbh - ctq relation needs to 
be converted to a Mbh — isph relation. One possible way is 
to estimate spheroid luminosities of the galaxies observed by 
T04 from their central velocity dispersions^ with the Faber- 
Jackson relation at 2 ~ 0.4 in the rest-frame B band (Ziegler 

^ According to T04, the velocity dispersions within an aperture 
used by T04 (uap) are converted to the central velocity dispersions 
(cro) as (TO = l.ltTap. 
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Figure 14. Upper panel: The Mbh — /'sph relations at z = 
(MH03; black line) and 0.35 (grey line) are shown. For the relation 
ai z = 0.35, passive evolution of Lsph is taken into account. The 
data points indicate the galaxies observed by T04 (see text for 
details). Lower panel: Shaded region indicates a spheroid-BHMF 
and its uncertainty obtained from the COMBO-17 LF at z = 
0.35 and the Mbh — /'sph relation at 2 = 0, which seems to 
be followed by the T04 galaxies. Solid lines show a BHMF at 
z = 0.25 obtained with a Mbh — ^sph relation considering a 
passive luminosity evolution between 2 = and 0.25. The upper 
and lower bounds are determined in the same way as those in 
FigureQSl The AGN-BHMF at 2 = 0.35 is overplotted with open 
circles. The error bars are calculated in the same way as those in 
Figure [T^ 



et al. 2005) and plot them against Mbh- The results are 
shown with solid circles in the upper panel of Figure 1141 
Although the sample size is small and there is a substantial 
scatter, the distribution of the data points suggests that the 
actual Mbh — isph relation followed by these spheroids lies 
closer to the local relationship than that allowing for a pas- 
sive evolution of Lsph between 2 = and 0.37. In the lower 
panel of Figure [Til a BHMF derived from the COMBO-17 
LF at z — 0.35 using the Mbh — isph relation at 2 = 0, which 
seems to be followed by the T04 galaxies, is indicated with 
shaded region. This BHMF is compared with that at a lower 
redshift (2 = 0.25) derived with a A/bh — isph relation con- 
sidering a passive luminosity evolution between 2 = and 
0.25. The latter BHMF is indicated with sohd lines. The un- 
certainties of these spheroid-BHMFs are indicated with the 
width of the shaded region or the separation of the lines; 
the upper and lower bounds are defined in the same way as 
those in FigureCSl The AGN-BHMF at 2 = 0.35 is overplot- 
ted with open circles. The error bars are calculated in the 
same way as those in Figure [T^ Although the uncertainties 
of the BHMFs are large, this comparison suggests that the 
BHMF obtained with the Mbh — Lsph relation followed by 
the T04 galaxies is slightly more biased to the massive end. 
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This trend could be seen more clearly when the BHMF is 
compared with that at 2: = 0. This apparently suggests a 
growth of SMBH towards higher redshift, which is unlikely 
in practice (see also Robertson et al. 2005). More galaxies 
need to be investigated to examine the correlation between 
Mbh and or L^ph at high redshift. 

4.3 Future perspective 

If much fainter end of early-type galaxy LF was determined 
from observations, BHMF could be probed further down 
to the low mass end [Mbh ^ 10* Mq). Since a redshift 
evolution of AGN-BHMF is suggested to be fast at inter- 
mediate redshifts in this mass range (Figure IIUI see also 
Merloni 2004), it would be interesting to look at the coun- 
terpart in spheroid-BHMFs in order to put much stronger 
constraints on the co-evolution of AGN and spheroid, par- 
ticularly down-sizing effects of their evolutions. Currently, 
galaxy LFs are not well constrained down to such low lu- 
minosities, especially at cosmological distances. Even at low 
redshifts, the faint end slope of an LF tends to be fixed to 
a certain value in fitting a Schechter function to the data 
and deriving L* and <j!>*. Much deeper data (e.g., ~ 2 mag 
deeper than the COMBO-17 limit) are essential to address 
the faint end of LF out to 2 ~ 1. Keeping a survey area 
similar to or even wider than COMBO-17 is the key to de- 
riving reliable LFs and to estimate cosmic variance. Multi- 
band photometry would be required to obtain photometric 
redshifts in good accuracy down to the faint end. In addi- 
tion, high spatial resolution images would enable one to di- 
rectly measure spheroid luminosities of galaxies. Although 
it is highly expensive to take data sets satisfying all these re- 
quirements, one promising candidate for this challenge is the 
COSMOS survey: a 2 square degree field is surveyed with 
the HST/ACS in the Jgi4 band down to 27.8 mag {5a) in 
the AB magnitude (cf. galaxies with ii 5C 24 mag are used to 
study LFs in COMBO-17) and also with Subaru/Suprime- 
Cam in the BVr'i'z' bands (Taniguchi et al. 2005). 
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